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Abstract 

We report a direct and robust calculation, free from ergodic problems, of the non- 
uniform-to-uniform (stripe) transition line of noncommutative by means of an exact 
Metropolis algorithm applied to the first non-trivial multitrace correction of this theory 
on the fuzzy sphere. In fact, we reconstruct the entire phase diagram including the Ising, 
matrix and stripe boundaries together with the triple point. We also report that the 
measured critical exponents of the Ising transition line agrees with the Onsager values in 
two dimensions. The triple point is identified as a termination point of the one-cut-to- 
two-cut transition line and is located at (6, c) = (—1.55,0.4) which compares favorably 
with previous Monte Carlo estimate. 

1 Introduction 

Noncommutative scalar phi-four theory is a two-parameter model which enjoys three 
stable phases: i) disordered (symmetric, one-cut, disk) phase, ii) uniform ordered (Ising, 
broken, asymmetric one-cut) phase and iii) non-uniform ordered (matrix, stripe, two- 
cut, annulus) phase. This picture is expected to hold for noncommutative/fuzzy phi-four 
theory in any dimension, and the three phases are all stable and are expected to meet 
at a triple point. The non-uniform ordered phase [27] is a full blown nonperturbative 
manifestation of the perturbative UV-IR mixing effect |26| which is due to the underlying 
highly non-local matrix degrees of freedom of the noncommutative scalar field. 

The phase structure in four dimensions was discussed using the Hartree-Fock approx¬ 
imation in [T] and studied by means of the Monte Carlo method, employing the fuzzy 
torus [3| as regulator, in [2|. 
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In two dimensions the theory is renormalizable [1]. The regularized theory on the 
fuzzy sphere reads 

S = Tr(a$[L„[L„$]]+6$2 + c^>4). (pi) 

The Laplacian A = Ca^a defines the underlying geometry, i.e. the metric, of the fuzzy 
sphere in the sense of UM- It is found that the collpased parameters are 

b = bN-^/^ja , c = cN-^ja'^. (1.2) 

The above phase structure was confirmed in two dimensions by means of Monte Carlo 
simulations on the fuzzy sphere in [9l[T0]. The phase diagram is shown on figures ([7]). 
Both figures were generated using the Metropolis algorithm on the fuzzy sphere. In the 
first figure coupling of the scalar field to a U{1) gauge field on the fuzzy sphere is 
included, and as a consequence, we can employ the U{N) gauge symmetry to reduce 
the scalar sector to only its eigenvalues. In the second figure an approximate Metropolis 
algorithm, i.e. it does not satisfy detailed balanced, is used. 

The problem of the phase structure of fuzzy scalar phi-four was also studied by means 
of the Monte Carlo method in |11II141[25] . The analytic derivation of the phase diagram of 
noncommutative phi-four on the fuzzy sphere was attempted in |15II2np24ll32] . The related 
problem of Monte Carlo simulation of noncommutative phi-four on the fuzzy torus, and 
the fuzzy disc was considered in m.m, and [22] respectively. For a recent study see [23] . 

In this paper, we are interested in studying by means of the Monte Carlo method the 
first non-trivial multitrace matrix model, quartic in the scalar field, which approximates 
noncommutative on the fuzzy sphere. The multitrace approach was initiated in [ISUIG] . 
See also [23] for a review and an extension of this method to the noncommutative Moyal- 
Weyl plane. For an earlier approach see |32j and for a similar more non-perturbative 
approach see [IM]. The multitrace expansion is the analogue of the Hopping parameter 
expansion on the lattice in the sense that we perform a small kinetic term expansion, i.e. 
expanding in the parameter a of m, while treating the potential exactly. This should 
be contrasted with the small interaction expansion of the usual perturbation theory. This 
technique is expected to capture the matrix transition between disordered and non-uniform 
ordered phases with arbitrarily increasing accuracy by including more and more terms in 
the expansion in a. From this we can then infer and/or estimate the position of the triple 
point. Capturing the Ising transition, and as a consequence the stripe transition, is more 
subtle and is only possible, in our opinion, if we include odd moments in the effective 
action and do not impose the symmetry <I> —)• — <I>. 

The effective action obtained in the multitrace approach is a multitrace matrix model, 
depending on various moments = TrM” of an A x A matrix M, which to the lowest 
non-trivial order is of the form 


V 


BTrM^ + CTrM‘^ + D 


2 


TrM^ 


+ B {TrMf+ C'TrMTrM^ + D {TrM)‘^ + XTrM‘^{TrMf + .... (1.3) 
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The parameters B and C are shifted values of b and c appearing in (11.11) . The primed 
parameters depend on a. The second line includes terms which depend on the odd mo¬ 
ments mi and m 3 . By diagonalization we obtain therefore the N eigenvalues of M as our 
independent set of dynamical degrees of freedom with an effective action of the form 

Seff = 

-b 


The coefficients U 2 ,i, and ^ 2,2 are given by the following two competing calculations 
found in [15] (Model I) and [21] (Model II): 

3 

U 2 ,i = -1 , U 44 = - , V 2,2 = 0 , Model I 

U 2,1 = -bl , U 44 = 0 , U 2,2 = ^ , Model II. (1.5) 

O 

The result found in [21] agrees with the non-perturbative result of m and the corrected 
result of m- This can also be confirmed by means of Monte Carlo. The first model in 
the commutative limit N — 00 is therefore a scalar theory on the sphere modulo 
multi-integral terms. In here, we will study both models by means of Monte Carlo and 
show that the first model, though incorrect, sustains the uniform ordered phase. The 
second model will sustain the uniform ordered phase only if we add to it higher order 
multitrace corrections. 

Since these models depend only on N independent eigenvalues their Monte Carlo 
sampling by means of the Metropolis algorithm does not suffer from any ergodic prob¬ 
lem and thus what we get in the simulations is really what should exist in the model 
non-perturbatively. This should be contrasted with the Monte Carlo simulation (via 
Metropolis, Hybrid Monte Carlo or other method) of (|l.ll) which suffers from severe er¬ 
godic problems which do not allow us easy and transparent access to the stripe transition 
and the triple point [nHans]. The model I, which sustains the uniform ordered phase, 
suffers, however, from critical slowing down, for values of N of the order of >60, and 
thus the use of the Wolf algorithm [M] would have been more appropriate. 

Some of our results in this article include: 

• The phase diagram of model I contains the three phases discussed above. The critical 
boundaries are determined and the triple point is located. 

• The uniform ordered phase exists in the model I only with the odd terms included. 
If we assume the symmetry M —> —M then the second line of (|1.3p becomes 
identically zero and the uniform ordered phase disappears. This is at least true in 
the domain studied in this article which includes the triple point of fuzzy on the 
fuzzy sphere and extends to all its phase diagram probed in pno]. 


^( 6 A 2 -b cA^) 


-A, 




-V2, 








24A^2 


V2,2 






(1.4) 
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• The delicate computation of the critical exponents of the Ising transition is discussed 
and our estimate of the critical exponents i/, a, 7 , /3 agrees very well with the Onsager 
values [ 33 ]. 

• The phase diagram of model II, with or without odd terms, does not contain the 
uniform ordered phase. 

• The one-cut-to-two-cut transition line does not extend to the origin, i.e. to C = 0, 
in the model II which gives us an estimation of the triple point in this case. 

• In the model II without the odd terms the termination point can be computed 
analytically from the requirement that the critical point I?* remains always negative 
and the obtained result {B,C) = (0,1/12) agrees with Monte Carlo. 

• In the model II with odd terms the termination point is found to be located at 
{B,C) = (—1.05,0.4). This is our measurement of the triple point. 

• In all cases the one-cut-to-two-cut matrix transition line agrees better with the dou¬ 
bletrace matrix theory than with the quartic matrix model. The doubletrace matrix 
theory is given hy D ^ 0 while all primed parameters are zero. 

• The model of Crosse-Wulkenhaar is also briefly discussed. 

This article is organized as follows: 

1. Section 2: The Multitrace Matrix Models. 

2. Section 3: Exact Solutions. 

• The Pure Real Quartic Matrix Model. 

• The Doubletrace Quartic Matrix Model. 

3. Section 4: Algorithm. 

4. Section 5: Monte Carlo Results. 

• General Remarks. 

• Monte Carlo Tests of Multitrace Approximations. 

• Phase Diagrams. 

• Gross-Wulkenhaar Model. 

5. Section 5: Conclusion. 

We also include appendices for the benefit of interested readers and to make the presen¬ 
tation as self-contained as possible. 


2 The Multitrace Matrix Models 


Our primary interest here is the theory of noncommutative on the fuzzy sphere 
given by the action 



( 2 . 1 ) 
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The Laplacian is A = [L^, [Lq, ...]]. Equivalently with the substitution <1> = M.ly/2'K6, 
where M. = Yl!ij=i ^ij\i >< j\i this action reads 


The parameters an 


S 


= Tr (^aMAM + bM^ + . 


1 



A 1 
4!^' 


In terms of the matrix M the action reads 


( 2 . 2 ) 


(2.3) 


S[M] = r‘^K[M]+ Tr + cM^]. 


(2.4) 


The kinetic matrix is given by 

- T+MTM - + EM‘^ 

The matrices T, Ts and E are given by 


K[M] = Tr 


(2.5) 


/ ITl 1 

(r3)zm = I5lm , (r)zm = W (m - 1)(1 - ^ ^ )^Zm-l , (-E)zm = (^ “ (2.6) 

The relationship between the parameters a and is given by 

= 2aN (2.7) 

We start from the path integral 


Z = 


J dM exp ( — S[Mf) 

j dAA2(A) exp(^-rr(6A2 + cA^)^ j dU exp . (2.8) 


The second line involves the diagonalization of the matrix M (more on this below). The 
calculation of the integral over U G U{N) is a very long calculation done in [151124] . The 
end result is a multi-trace effective potential given by 


Seff 


+ 


+ cXf) - ^ X] 

i 

'jV2,i ^(Ai - Xjf + ^U4,1 ^(Ai - XjY 
^ i¥=j 


24A2 


V2,2 


+ ••• • 

i¥=j 

(2.9) 


^The noncommutativity parameter on the fuzzy sphere is related to the radius of the sphere by 9 

2AVVAf2 _ p 
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The coefficients v will be given below. 

We will assume now that the parameters b, c and scale as 


a = 


, h = 


c = 


r = 


JSf&a. ’ ]\[Sb ’ Ar<5c ’ ]\[&r 

Further, we will assume a scaling 5x of the eigenvalues A, viz 




( 2 . 10 ) 


( 2 . 11 ) 


It is easy to convince ourselves that in order for the above effective potential to come out 
of order N'^ we must have the following values 


5, = -1 - 25a , = 1 - 2,5a , 5c = 1 - 45 a , = -25a. (2.12) 

From the Monte Carlo results of [aio], we know that the scaling behavior of the param¬ 
eters h and c appearing in the above action on the fuzzy sphere is given by 

5fe = ^ , 5c = 2. (2.13) 

By substitution we obtain the other scalings 

h = -\, Sa = -^, 6r = (2.14) 

The problem (j2.9j) is a generalization of the quartic Hermitian matrix potential model. 
Indeed, by dropping odd moments, this effective potential corresponds to the matrix model 
given by 


C = Co + ACo- 

The classical potential and the even correction AVq are given by 

Co = bTrM'^ -h cTrM"^. 


(2.15) 


(2.16) 


ACo = F'TrM"^ -h e'TtM^ + D 


TrM^ 


The coefficients F , E and D are given by 


/ aiV^U2,i a^N'^Vi^i 2r]a?N‘^ 

r = - , hj = - , U = -. 

2 ’ 6 ’ 3 

The strength of the multi-trace term r] is given by 


r] = V2,2 - ^^4,1. 


(2.17) 


(2.18) 


(2.19) 


By including terms which involve the odd moments we get the effective potential 


C 


Co + ACo + AC. 
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( 2 . 20 ) 












( 2 . 21 ) 


The extra contribution and its coefficients are given by 

AF = b'{T rMf+ C'TrMTrM^ + d'{T rMf + XTrM‘^{TrMf. 


B = 


aN 

—V2,l 


C = 


2a‘^N^ 

— - —r;4,i , D = 


2a^ , 



( 2 . 22 ) 


The coefficients r; 2 ,i, ^ 4,1 and ^ 2,2 are given by the following two competing calculations 
found in [15] (Model I) and jUj (Model II): 

3 

U 2 ,i = -1 , U 4 ,i = - , U 2,2 = 0 , Model I 

V2,1 = +1 , V4:,i = 0 , V2,2 = ^ , Model II. (2.23) 

O 

The difference in the sign of U 2 ,i is probably a typo on the part of m while the discrepancy 
in the values of U 4 _i and ^ 2,2 is more serious and is discussed in [Mj- The result found 
in [231 agrees with the result of m given by their equation (2.39). The work [16] contains 
the correct calculation which agrees with both the results of m and [23] . 

The one-cut-to-two-cut transition line in the model (|2.15l) is given by the exact result 

on 


= -I-.,. -2jc+ . (2.24) 

2 '' « 3y'c + f t. 4.1 

As point out in [23] this result is new. For a generalization of this result see |17j . This 
critical value 6 * is negative for 

c>-(477-U4 ,i). (2.25) 
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3 Exact Solutions 

In this section we will give a brief description of the exact solutions of the real quartic 
matrix model BTrM"^ + CTrM^ and the doubletrace real quartic matrix model BTrM‘^ + 
CTrM"^ + D{TrM‘^f. 


3.1 The Pure Real Quartic Matrix Model 


The phase structure of the pure real quartic matrix model is studied for example 
in |28H31j . In here we will summarize some of the salient results. 

The basic model is given by 


V = BTrM^ + CTrM"^ 

9 4 ' 


(3.1) 


B 


a r-E 

9 ’ 4ff- 


(3.2) 


There are two phases in this case: 
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Disordered phase (one-cut) for g > Qc'. 


p{X) = -^{2CX^ + B + CS^)VS‘^ - A2 

Ntt 

= —(-A^ — 1 + r^)-\/4r2 — A2. 

gn 2 

(3.3) 

-2r <X<2r. 

(3.4) 

1. 

r = — 0 . 

2 

(3.5) 

= -^i-B + ^yB^ + UNC) 

= ^(1 + Vi + 

(3.6) 

Non-uniform ordered phase (two-cut) for g < g^. 

P(\) = ^\/(A2-5f)(«|-A2) 

- 2^1 

(3.7) 

r- < A < r+. 

(3.8) 

r_ = (5i , r+ = 62 . 

(3.9) 

= ^{-B^ 2 VnC) 

= 2(lT\/ff). 

(3.10) 


Critical point: A third order transition between the above two phases occurs at the 
critical point 


Qc = I ^ BI = ANC ^ Bc = -2VNC. 


(3.11) 


Specific heat: The behavior of the specific heat across the matrix transition provides 
also a powerful result against which we can calibrate our algorithms and Monte Carlo 
simulations. In terms of B = B/Be the specific heat reads in the two phases of the theory 
as follows 


a 

C, 


= B = B/Be < -1 

1 2B^ B -rf - 

= - +-(2^2-3 VS2 +3 

4 27 27 ^ ^ 


B > -1. 


(3.12) 



Uniform ordered phase: The real quartic matrix model admits also a solution with 
TrM 7 ^ 0 corresponding to a possible uniform-ordered (Ising) phase. This like 

solution can appear only for negative values of the mass parameter and it is constructed, 
for example, in [28] , It is, however, thought that this solution can not yield to a stable 
phase without the addition of the kinetic term to the real quartic matrix model. 

The density of eigenvalues in this case is given by 

p{z) = -^—{2Cz‘^ -|- 2aCz + B + 2Ca‘^ + (3.13) 

ttN 

This is a one-cut solution centered around r in the interval [cr — r, cr -|- r] where a and r 
are given by 

+ 2^/B^ - 15NC) , = r^{-2B - 2yJB'^ - IbNC). (3.14) 

iUC/ loO 

This solution makes sense only for 

B <Bc = -^/T5^/NC. (3.15) 

3.2 The Doubletrace Quartic Matrix Model 

The doubletrace real quartic matrix model is given by the multitrace matrix model 
(m with all odd moments set to zero, viz 


V = BTrM^ + CTrM"^ + D(TrM^)^. (3.16) 

The scaling of the parameters is given by 

B = , C = CN-‘^ , D = DN-^. (3.17) 


The phase structure of this model is very similar to the phase structure of the pure 
real quartic matrix model outlined in the previous section. See for example [2Tj- The two 
stable phases are still given by the disordered (one-cut) phase and the non-uniform-ordered 
(two-cut) phase separated by a deformation of the line .B* = —2^/c given by 


B, = -2Vd - (3.18) 

For a generalization of this result see m- 

Another important result for us here is the existence of a termination point in the 
model of [M] since the critical line does not extend to zero. Indeed, in order for the 
critical value B^, to be negative one must have C in the range 


c>a = 


2r/d^ 


a 

12 ' 


Thus the termination point is located at (for a = 1) 

(B,C) = (0,1/12). 


(3.19) 


(3.20) 
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4 Algorithm 


We start from the potential and the partition function 


V = Tr {BM‘^ + CM"^) + D {TrM'^) ^ 

+ b' {TrMf + C'TrMTrM^ + d'{T rM)^ + A'TrM^{TrMf 


(4.1) 


Z = j dM exp(-y). (4.2) 

The relationship between the two sets of parameters {a, b, c} and {B, C, D} is given by 

(4.3) 


aN‘^V2,i „ 2r]a^N‘^ 

B = b-\ --—^ , C = c H--—^ , D = ---. 


6 


The collpased parameters are 


5 B ~ dv 2 ,i ^ C d‘^VA,i ^ 2r]a?N 

B = — = b+—- , C = ^ = c + —— , D = --—. 

jV2 2 A'^ 6 3 


(4.4) 


Only two of these three parameters are independent. For consistency of the large N 
limit, we must choose d to be any fixed number. We then choose for simplicity a = 1 or 
equivalently D = — 2r/A^/li^. The other parameters are 


aN 


B = «,i , C = - ^ 


2a2Ar2 , 2a2 , do^A^ 

U4,l , D = - ^'^2,2 , A = - - - V2,2- 


We can now diagonalize the scalar matrix M as 


M = UAU-\ 


We compute 


6M = u(^6A + [U-^6U,A]^U-\ 


(4.5) 

(4.6) 

(4.7) 


Thus (with U ^5U = i6V being an element of the Lie algebra of SU(N)) 
Tr{5Mf = Tr{5Af + Tr[U-^5U,Af 


(4.8) 


i¥^i 


We count A^^ real degrees of freedom as there should be. The measure is therefore given 
by 


dM = dXi dVijdV*jy/det{metric) 

i i^j 


X{dX,\{dV,,dV* 

i y i^j 


(4.9) 


^The authors of mm chose instead a = 1. This should not make any difference to the Monte Carlo 
simulations. 


10 













We write this as 


dM = dAdUA‘^{A). 


(4.10) 


The dU is the usual Haar measure over the group SU(N) which is normalized such that 
f dU = 1, whereas the Jacobian A^(A) is precisely the so-called Vandermonde determinant 
defined by 


A2(A) = n(A,-A,)2. 

i>j 

The partition function becomes 

Z = j dA A2(A) exp ^ - Tr{BA^ + CA^) - D^TrA^^ ^ . 


We are therefore dealing with an effective potential given by 

i=l i=l i=l '' ^ ijij 


(4.11) 


(4.12) 


(4.13) 


We will use the Metropolis algorithm to study this model. Under the change A* —)■ \i + h 
of the eigenvalue A^ the above effective potential changes as Ueff —^ where 


AVi^h = BAS2 + CAS4 + D(2S2AS2 + ASi) + ASv^nd 

+ b'AS 2 + C'AS'^ + D'f {AS'^f + A 52 ) + A' I (Si + h)AS2 + hS2 


(4.14) 


The monomials Sn are defined by Sn = Yli ^i' while the variations ASn and A5vand are 
given by 


AS2 = h^ + 2hXi. 

(4.15) 

ASi = e/r^Af + AhXf + Ah^Xi + h\ 

(4.16) 

^5'vand — 2 ^ In 1 J- . 

jV* * ^ 

(4.17) 

AS2 = h^ + 2 hSi. 

(4.18) 

AS4 = (Si + h){ 3 hX‘f + 3 h^Xi + h^) + hS^. 

(4.19) 
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Monte Carlo Results 


5.1 General Remarks 

1. We use the statistics 2^ + 2^ x 2^ with P = 15 — 20 and P' = 5 with = 10 — 60 
and with the Jackknife method to estimate the error bars. We can even go further 
to = 100 and beyond but noticed that critical slowing down became a serious 
obstacle especially in the measurement of critical exponents. 

2. Our hrst test for the validity of our simulations is to look at the Schwinger-Dyson 
identity given for the full multitrace model (j 2 . 20 ll by 

< {2bTrM‘^ + AcTrM^ + 2 V 2 + 414) >= (5.1) 

The quartic and quadratic pieces V 2 and V 4 are such that 

AI/q + AF = F 2 + V4. (5.2) 

In other words, 

V 2 = f'TtM^ + b' {TrMf. (5.3) 


V 4 = e'TvM^ + D 


2 




+ C'TrMTrM^ + d'{T rMf + A'TrM^{TrMf .{5A) 


3. The second powerful test is to look at the conventional quartic matrix model with a = 
0, viz V = Vq. The eigenvalues distributions in the two stable phases (disorder(one- 
cut) and non-uniform order (two-cut)) as well as the demarcation of their boundary 
are well known analytically given by the formulas (13.31) . (|3.7p and (13.lip . 

4. Even the quartic multitrace approximation itself can be verihed directly in Monte 
Carlo in order to resolve the ambiguity in the coefficients v between [23] and US). 
We must have as identity the two equations 

<aj dUTr[La, UAp-^f >Vo=< -V 2 {A) >v, ■ (5.5) 


<\{a j dUTr[La, UAU-^f^ >vb=< -f4(A) + ^^^^(A) >y„ . (5.6) 

The coefficients v appear in the potentials V 2 and V 4 . The expectation values are 
computed with respect to the conventional quartic matrix model Vq = Vo(A). 

This test clearly requires the computation of the kinetic term and its square which 
means in particular that we need to numerically perform the integral over U in 
the term f dPTr[La, l7AU~^f which is not obvious how to do in any direct way. 
Equivalently, we can undo the diagonalization in the terms involving the kinetic term 
to obtain instead the equations 

< aTr[La, Mf >vb=< —V2 >Vo ■ 
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(5.7) 








(5.8) 



aTr[La,M\ 


>Vo=< -V4 + -V2 >Vo ■ 


Now the expectation values in the left hand side must be computed with respect to 
the conventional quartic matrix model Vq = Vb(M) with the full matrix M = UAU~^ 
instead of the eigenvalues matrix A. The expectation values in the right hand side 
can be computed either ways. 


In other words, the eigenvalues Metropolis algorithm employed in this article to 
compute terms such as < — V2 >vb and < — V4 + V 2/2 >Vo can not be used to 
compute the terms < aTr[La, >Vq and < (aTr[La, M]^)^/2 >Vo- We use instead 
the hybrid Monte Carlo algorithm to compute these terms as well as the terms 
< —V 2 >Vo and < — V 4 + V 2/2 >Vo in order to verify the above equations. This 
also should be viewed as a counter check for the hybrid Monte Carlo algorithiri^ 
since we can compare the values of < —V 2 >vb and < — V4 + V 2 /2 >v"q obtained 
using the hybrid Monte Carlo with those obtained using our eigenvalues Metropolis 
algorithm. We note, in passing, that the Metropolis algorithm employed for the 
eigenvalues problem here is far more efficient than the hybrid Monte Carlo applied 
to the same problem without diagonalization. But this we can obviously tolerate for 
testing purposes. 


5. The most detailed order parameter at our disposal is the eigenvalues distribution of 
the field/matrix M which behaves in distinct ways in various phases. This behavior 
mimics their behavior in the conventional quartic matrix model Vq, viz 

• The disorder (one-cut) phase is characterized by a single-cut eigenvalues distri¬ 
bution symmetric around 0 since in this phase < M >= 0. 

• The non-uniform order (two-cut) phase is characterized by an eigenvalues dis¬ 
tribution symmetric around 0 but with two disjoint supports since < M >= 
\J—hl2c 7 where 7 is any A^—dimensional idempotent, i.e. 7 ^ = 1. This appears 
for large values of c. 

• The uniform order (asymmetric one-cut) phase is characterized by a single¬ 
cut eigenvalues distributions centered around a non-zero value since < M >= 
y/—bj2c 1. This appears for small values of c. 

6 . The sepcihc heat defined with the respect to the multitrace potential is given by 


Cy=<V^>-<V>^. 


(5.9) 


The relation of this powerful and most difficult to measure second moment with the 
specific heat of noncommutative on the fuzzy sphere is discussed in the appendix. 

In any case, this specific heat is expected to approach the specific heat of the original 
noncommutative for large values of c, and as a consequence it can be used to 
locate the boundary between one-cut and two-cut as in the conventional quartic 
matrix model with a = 0 . 

^Or a counter check for the eigenvalues Metropolis algorithm depending on which algorithm is more trust¬ 
worthy. However, we firmly believe that the eigenvalues Metropolis algorithm used here is more robust on all 
accounts. 
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7. The magnetization and susceptibility are defined by 

m =< |TrM| > , y =< |TrMp > — < |rrM| . (5.10) 

The magnetic susceptibility will exhibit peaks in the second order phase transitions 
between disorder (one-cut) and uniform order (Ising) and between non-uniform order 
(two-cut) and uniform order (Ising). 

8. The total power and power in the zero mode are defined by 

Pt =< ^TrM^ > , Po=< {^TrMf > . (5.11) 

In the Ising (uniform order) phase we will have in particular the very distinguished 
signal P = Pq. 

5.2 Monte Carlo Tests of Multitrace Approximations 

It is quite obvious that resolving the ambiguity between the calculations of [l5] and [21] , 
summarized in equation (|2.23p . is straightforward in Monte Carlo. We only need to 
show that the two equations (|5.7p and (15.8p hold as identities in the correct calculation. 
However, this requires a different algorithm than the eigenvalues Metropolis algorithm 
used here. Indeed, to solve this problem we need to Monte Carlo sample, both the 
eigenvalues and the angles of the matrix M using the Metropolis or the hybrid Monte 
Carlo, the quartic matrix model 


Vo = bTrM‘^ + cTrM^. (5.12) 

Monte Carlo simulations of this model can also be compared to the exact solution outlined 
in section 3.1 so calibration in this case is easy. The detail of this simple exercise is 
reported in [24]. There, it is decisively shown that the calculation of [24] gives the correct 
approximation of noncommutative scalar <I>| on the fuzzy sphere. 

5.3 Phase Diagrams 

1. Model I: Model of |15| : 

• Ising: Some of the results for the Ising transition for this model are shown on 
table ©• The critical point is taken at the peak of the susceptibility. The be¬ 
havior of various observables is shown on figure m- The fit for the extrapolated 
critical value is given by 

C = 0.291(0).(-H)-h 0.104(1). (5.13) 

This transition can be confirmed to be between disordered and uniform-ordered 
by looking at the eigenvalues distribution. In the disordered phase we have 
one-cut symmetric around zero whereas in the uniform-ordered we have one-cut 
symmetric around H/2C. See figure ([2]). 
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• Critical Exponents: The calculation of the critical exponents of the above 
Ising transition is a very delicate exercise in Monte Carlo due to the known 
problem of critical slowing down and as a consequence the use of a different 
algorithm, for large values of N, such as the Wolf algorithm [M] is essential. 
For values of N less than N = 60 the current algorithm is sufficient. In any 
case, this lengthy calculation is reported elsewhere. Suffice to say here that 
the critical exponents obtained are consistent, within the best statistical errors, 
with the Onsager solution of the Ising model in two dimensions given by the 
celebrated values [33] 


u = l , /3 = l/8 , 7 = 7/4, a = 0 , r] = 1/4. (5.14) 

This has always been known to be true but this is the first Monte Carlo direct 
calculation of these critical exponents. 

• Matrix: Some of the results for the matrix transition between disorder and 
non-uniform order for this model are shown on table (|2|). The critical point 
is determined at the point where the eigenvalue distributions go from one-cut 
in the disorder phase to two-cut in the non-uniform phase. The splitting of 
the distribution is considered to have been occurred when the hight of the 
distribution at A = 0 is less than some tolerance Tol. We take Tol = 0.001. 
The behavior of the specific heat across this transition is effectively that of the 
pure quartic matrix model a = 0. A sample of the corresponding specific heats 
and eigenvalue distributions is shown on figure (|3|). The fit for the extrapolated 
critical value is given by 

C = 2.206(67).(-S) - 7.039(301). (5.15) 

• Stripe : This transition is quite difficult to observe in Monte Carlo even in this 
simplified setting which involves the sampling of N eigenvalues. We can observe 
this transition for medium values of C immediately above, but not too close to, 
the triple point. The transition point is taken at the point where we observe a 
jump or a discontinuity in the zero power Pq and the specific heat as seen on 
figure ([I]). 

Alternatively, we can approach the critical boundary by fixing the value of 
B and changing C starting from small values, i.e. inside the uniform ordered 
phase, until the curves for the total and zero powers start to diverge marking the 
transition to the non-uniform ordered phase. The signal we obtain in this way is 
quite clear and unambiguous as shown on figure (|5|) and some measurements are 
included in table (|3|). Since this is a very delicate transition we do not perform 
any extrapolation of the critical point and the critical boundary is given by the 
fit of the largest value of N. In any case we observe no strong dependence on 
N of the measured critical value C as seen on table ([3|). The stripe critical line 
is then approximated by the fit for N = 50 given by 

C = 0.154(22).(-S) -F 0.530(131) , = 50. (5.16) 
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• Triple Point and Phase Diagram: The location of the triple point is obtained 
from the intersection point of the Ising and matrix lines (j5.13p and (j5.15p re¬ 
spectively. Indeed, the measurement of these two lines is more robust than the 
measurement of the stripe line (15.161) . We get then immediately 

(5,(7) = (-3.73,1.19). (5.17) 

The phase diagram of the multitrace model of m is shown on figure (l6|). The 
Ising and matrix transition data points are not shown explicitly but we only 
include their extrapolated fits whereas the N = 50, 36 and 25 stripe data points 
are indicated explicitly. We observe that the matrix boundary is closer to the 
doubletrace theory than it is to the quartic matrix model. The stripe critical 
boundary is of course expected to be closer to the = 50 measurement. 

• Even Model: This is the model in which we set all odd moments to zero in 
the action. We get then the doubletrace model 

V = BTrM^+ CTrM‘^ + D(TrM^f. (5.18) 

The most fundamental property of this model, observed in Monte Carlo simula¬ 
tion, is the absence of the uniform ordered phase. Indeed, only the disorder and 
the non-uniform order phases exist in the phase diagram. The critical boundary 
is very close to the doubletrace critical line shown on figure Q which consists 
of two branches. Some precise measurements for the first branch are included in 
table ([4]). The turning point, towards the second branch, occurs around C ~ 0.1 
where the critical point —B^ becomes increasing, instead of decreasing, as we 
decrease C. 

2. Model II: Model of [24j : This model as pointed out previously is the correct 
approximation of noncommutative scalar <1>2 on the fuzzy sphere. However, this 
model is characterized by the absence of the uniform ordered phase and only the 
matrix transition line separating disordered and non-uniform ordered phases exists 
in the phase diagram. This fundamental result holds with and without odd terms. 
The role of the odd terms seems to be negligible and the two cases with and without 
odd terms are close. The doubletrace theory is also a very good approximation. A 
phase diagram is attached on figure (j6|). 

The second fundamental observation in this case is the existence of a termination 
point. The matrix critical line does not extend to the origin and terminates at 
a point around C = 0.083 in the case without odd terms, which agrees with the 
doubletrace theory prediction (|3.20l) . and at a point around C = 0.4 in the case with 
odd terms. This termination point is exhibited in Monte Carlo by the failure of the 
Schwinger-Dyson identity (|5.ip . 

In particular, we observe for N = bQ that C = 0.4 is the smallest value at which the 
disordered (one-cut) and non-uniform ordered (two-cut) phases are well defined. For 
C = 0.2 — 0.3 the non-uniform ordered phase can not be clearly observed whereas 
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for C < 0.1 both the disordered and the non-uniform ordered phases become in¬ 
discernible. It is therefore natural to identify the triple point with the termination 
point C = 0.4. Our estimation of the termination point is given by 

(.B,C) = (-1.05,0.4). (5.19) 


5.4 Grosse-Wulkenhaar Model 


The multitrace approach can also be applied to a regularized noncommutative 4>2 on 
the Moyal-Weyl plane in the matrix basis [24j with action given by 


S = Tr 


N 


+ a(EM^ + ^/UT+MTM 




(5.20) 


Two cases are of importance to us here: 


1. The noncommutative theory without a harmonic oscillator term. In this case the 
effective action takes the form 

Seffe = bTrNM^+ cTrNM'^ + d{TrNM^f + bi{TrNMf+ ci{TrNM)'^ 

+ diTri^M‘^{Tr]^M)‘^ + eTrj^MTrj^M^. 


(5.21) 


The parameters are given by 

, aN 

b= -^- 

2 2 


u a^N , a 
c = — - , d = — 


N 24 


2 

12 


a a 

^ “ “2 ’ 24iV2 


= — 


e = 


(5.22) 


If we assume the symmetry M 
we end up with the action 


12 Ar ' ' 6 

—M then all odd moments vanish identically and 


2^2 


S’effe = bTrNM^ + cTrNM^ + d{TrNM^) 


(5.23) 


2. At the self-dual point we have = 1, and thus ^/uj = 0, and as a consequence the 
effective action reduces to the multitrace model 


S’effe = bTrNM^ + cTr^M^ + d{TrNM^f. 


The parameters b, c and d are given by 



c 


u 

N 


a?N 

~W 


d 



(5.24) 


(5.25) 


Both the actions (I5.23|) and (j5.24p do not contain odd moments and thus the corresponding 
phase diagrams are expected to not contain the uniform ordered phase with all matrix-like 
behavior as consequence. 
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c 

AT = 10 

N = 25 

N = 36 

iV = 50 

B extrapolated 

0.3 

-0.71 

-0.69 

-0.68 

-0.68 

-0.672(2) 

0.5 

-1.44 

-1.39 

-1.38 

-1.38 

- 1.361(3) 

1.0 

-3.21 

-3.13 

-3.11 

-3.10 

-3.073(2) 

1.2 

-3.9 

-3.82 

-3.8 

-3.79 

-3.763(2) 


Table 1: The Ising transition points for iV = 10 — 50. These are determined at the peak of the 
susceptibility (discontinuity in the specihc heat). The search step is 0.01. 


C 

iV = 10 

N = 25 

N = 36 

iV = 50 

B extrapolated 

2.0 

-4.925 

-4.475 

-4.325 

-4.225 

-4.090(37) 

2.5 

-5.225 

-4.725 

-4.575 

-4.475 

-4.321(31) 

3.0 

-5.525 

-4.975 

-4.875 

-4.725 

-4.576(36) 

4.0 

-6.025 

-5.525 

-5.225 

-5.175 

-4.993(83) 


Table 2: The matrix transition points for iV = 10 — 50. These are determined at the point 
where the eigenvalue distribution splits which is taken at the value of B where the distribution 
drops below 0.001 at zero. The search step is 0.025. 


B 

N 

= 

10 

N 

= 

25 

N 

= 

36 

N 

= 

50 

-9.0 

1.95 

T 

0.55 

1.8 

T 

0.4 

1.75 

T 

0.15 

1.95 

T 

0.15 

-8.0 

1.85 

T 

0.15 

1.65 

T 

0.25 

1.8 

T 

0.2 

1.8 

T 

0.2 

-7.0 

1.6 

T 

0.4 

1.45 

T 

0.05 

1.55 

T 

0.15 

1.6 

T 

0.3 

-6.5 

1.65 

T 

0.15 

1.45 

T 

0.05 

1.45 

T 

0.05 

1.45 

T 

0.15 

-6.0 

1.65 

T 

0.15 

1.35 

T 

0.05 

1.55 

T 

0.05 

1.5 

T 

0.2 

-5.5 

1.55 

T 

0.05 

1.35 

T 

0.05 

1.35 

T 

0.15 

1.35 

T 

0.05 

-5.0 

1.45 

T 

0.05 

1.45 

T 

0.05 

1.25 

T 

0.15 

1.4 

T 

0.1 


Table 3: The non-uniform-to-uniform transition points C for N = 10—50. These are determined 
at the discontinuity of the zero power. 
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c 

N=10 

iV = 17 

N = 25 

N = 36 

iV = 50 

B 

0.5 

-4.325 

-4.225 

-4.025 

-3.875 

-3.825 

-3.738(72) 

1 

-4.375 

-4.175 

-4.025 

-3.875 

-3.775 

-3.685(54) 

2 

-4.925 

-4.675 

-4.475 

-4.325 

-4.225 

-4.098(50) 

3 

-5.475 

-5.275 

-5.025 

-4.825 

-4.725 

-4.601(84) 

5 

-6.525 

-6.275 

-6.025 

-5.725 

-5.625 

-5.479(108) 


Table 4: The matrix transition points for N = 10 —50 in model I without odd terms. In this case 
only this transition exists and extends to a turning point in accordance with the doubletrace 
theory. The search step is 0.025. 

6 Conclusion 

A Monte Carlo study of the multitrace quartic matrix model of m , which is claimed to 
be the first non-trivial correction to noncommutative on the fuzzy sphere, is presented. 

This model does not suffer from the severe ergodic problems encountered in the simulations 
of noncommutative on the fuzzy sphere and the Metropolis algorithm is very effective 
in probing the entire phase space. In particular, Monte Carlo measurement of the one- 
cut-to-two-cut and the Ising transition lines as well as a direct Monte Carlo measurement 
of the non-uniform-to-uniform transition line are performed. The odd terms in the action 
which are dropped in m do play the central role in generating the Ising phase and the 
non-uniform-to-uniform transition line and thus a triple point. A quantitative sketch of 
the phase diagram and the triple point is outline. 

The closely related multitrace quartic matrix model of [24j . which is the correct ap¬ 
proximation of noncommutative scalar 4>2 on the fuzzy sphere, is also considered in this 
article where it is shown that the one-cut-to-two-cut transition line does not extend to 
the origin and terminates at a point consistent with the triple point of noncommutative 

on the fuzzy sphere [9]. 

We also commented in this article on the Crosse-Wulkenhaar model which is examined 
using a combination of the multitrace technique and the Monte Carlo method. At this 
order of the multitrace approximation the two models obtained in this case do not exhibit 
the Ising phase and the non-uniform-to-uniform transition line. 

A Susceptibility and Specific Heat 

A.l Susceptibility 

We consider on the fuzzy sphere coupled to a constant magnetic field H given by 
the action 


5 = Tr{a^La, [La, $]] + + H<^). 


(A.l) 
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The magnetization and the susceptibility are defined by 


magnetization = 


— < Tr^ > 
N 


' NdH 


In Z. 


susceptibility = < (Tr<h)^ > — < Tr$ 


dm 


\nZ 
d 


= —N— —magnetization. 

dH 


On the fuzzy sphere we have 


2^ r r 1 a rr. ^ 

Xa = -J^La , [Xa,Xb\ = —€abcXc , 0 = — , Tr = 

The regularized noncommutative plane is then defined by 


(fx. 


1 If 

X3 = R , [xi,X 2 ] =i9 , di = = -Q^ijXj , / (fx = 2-K6Tr. 

We have ei 2 = 1- The above action becomes, including a rescaling of the field 
YTVa/^^h, given by the equation 

S = 2TT6Tr (^-(j)didi(j) + -rrf(jf + + h(j)). 


m^ = — A = h = 

11 ^ , /V Ar_9 7->9 ’ 


N H 


aR^ ’ Na?R^ ' V 27ra 2i?2’ 

The commutative limit is 6 —)■ 0. By using a lattice in this limit we have 

S = + 2>^(Pn + h(pn) ■ 


We compute in this limit on the lattice 


magnetization = ^ Trd* > 




47ri?2 A/2 


< 


\r2 Y1 


> . 


The volume of the lattice must be equal to the area of the sphere, viz = AttR‘ 
we compute 


susceptibility = < (Tr<f>)^ > — < Tr<h 




< 




> - < 


>')• 


(A.2) 

(A.3) 

(A.4) 

(A.5) 
-^<p = 

(A.6) 

(A.7) 

(A.8) 

(A.9) 
*. Also 

(A.IO) 
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A.2 Specific Heat 

The specific heat is defined by 

= <S^ > - <S>^ . (A.ll) 

The inverse temperature is introduced in the usual way as 

Z = y (iMexp(—/35[M]). 

The calculation of the effective potential proceeds as before with the replacement a 
The partition function in the quartic multitrace approximation is 


(A.12) 
—a/3. 


Z = J dA A\A) exp {-m)+Pj dA A^A) exp ( - pVo) ( - V 2 ) 

+ j dA A2(A) exp ( - pVo) ( - 14 + Ivi) . (A.13) 

A straightforward calculation yields 

Cy = <{V + > - <{V + V4) >2 -2 < (14 + 2Vi + 2 V 2 V 4 ) > . (A.14) 

The last term could make this approximation of the specific heat negative. This actually 
happens in the approximation of |15] . 

Acknowledgment: This research was supported by CNEPRU: ’’The National (Al¬ 
gerian) Commission for the Evaluation of University Research Projects” under contract 
number D01120130009. 
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Figure 1: Some observables of the multitrace model of na across the disorder-to-uniform 
transition. 
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Figure 2: The eigenvalues distribution of the matrix M in the multitrace model of [I5] with 
C* = 1.0 across the disorder-to-uniform transition. 
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Figure 3: The specific heat and the eigenvalues distribution of the matrix M in the multitrace 
model of [15] across the disorder-to-non-uniform transition. 
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Figure 4: The total and zero powers and the specihc heat, as functions of B, of the multitrace 
model of [U] across the non-uniform-to-uniform transition. 
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Figure 5: The total and zero powers, as functions of C, and the eigenvalues distribution of the 
multitrace model of [15] across the non-uniform-to-uniform transition. 
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Figure 6: The phase diagrams of the multitrace models of jTS] and [24]. Model I; The Ising 
and matrix transition data points are not shown but only their extrapolated fits are included 
whereas the N = 25, N = 36 and = 50 stripe data points are indicated explicitly. Model II: 
the triple point is identified as the termination point located at {B,C) = (—1.05,0.4) which is 
to be compared with the doubletrace prediction at (0,1/12). 
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Figure 7: The phase diagram of noncommutative phi-four theory on the fuzzy sphere. In the 
first figure the fits are reproduced from actual Monte Carlo data izg. Second figure reproduced 
from [9] with the gracious permission of D. O’Connor. 
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